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The ability to measure and manipulate complex particles in the solid state is a cornerstone of 
modern condensed-matter physics. Typical excitations of a sea of electrons, called quasiparticles, 
have properties similar to those of free electrons. However, in recent years exotic excitations with 
very different properties have been created in designer quantum materials, including Dirac fermions 
in graphene^ and fractionally-charged quasiparticles in fractional quantum Hall systems^. Here we 
report signatures of a new quasiparticle - the Majorana fermion in .Josephson junctions consist- 
ing of two superconducting leads coupled through a three-dimensional topological insulator'^ We 
observe two striking departures from the common transport properties of Josephson junctions: a 
characteristic energy that scales inversely with the width of the junction, and a low characteristic 
magnetic field for suppressing supercurrent. To explain these effects, we propose a phenomenologi- 
cal model in which a one-dimensional wire of Majorana fermions is present along the width of the 
junction, similar to a theoretical prediction by Fu and Kane*. These results present an opening 
into the investigation of Majorana fermions in the solid state and their exotic properties, including 
non-Abelian statistics^, a suggested basis for fault-tolerant quantum computation^. 



The Majorana fermion, a charge-neutral particle that 

is its own antiparticlc, was proposed theoretically almost 
75 years ago^. In the original context of high-energy 
physics, neutrinos have been proposed to be Majorana 
fermions, a designation that could help explain the dom- 
ination of matter over antimatter in our universe^ . Sepa- 
rately, electronic excitations in certain condensed matter 
systems have recently been predicted to act as Majorana 
fermions'^. One such system is a three-dimensional topo- 
logical insulator (TI) where superconducting correlations 
between particles are introduced, producing a "topolog- 
ical superconductor" with properties similar to a p + ip 
superconductor^. Such superconducting correlations can 
be induced by doping a TI with Cu, where deviations 
from conventional superconductivity have been seen^. 
Another way to create a topological superconductor is 
via the proximity effect, wherein a conventional super- 
conductor is placed in contact with a TI^^. When two 
superconductors are connected by a TI, the TI "weak 
link" supercondiicts due to its proximity to the super- 
conducting leads. This produces a Josephson junction 
(JJ) but with several important distinctions compared 
to a conventional J J, where the weak link is typically an 
ordinary metal or insulator. Fu and Kane have predicted 
[4] a one-dimensional (ID) mode of Majorana fermions 
at the interface between a conventional superconductor 
and a superconducting topological surface state. Hence, 
JJs formed with a TI weak link are expected to have two 
ID modes at the two superconductor-TI interfaces (ar- 
rows in Fig. la), which fuse to form a ID wire of Majo- 



rana fermions (shown in purple in Fig. la) running along 

the width of the device^. The energy spectrum of these 
Majorana fermions is characterized by states within the 
superconducting gap, which cross at zero energy when 
the phase difference (p between the two superconducting 
leads is n (solid purple line. Fig. lb). 

To probe this exotic state, recent experiments have in- 
vestigated transport in TI JJs, finding good agreement 
with conventional J J behavior^^"!^. Two characteris- 
tic properties are typically reported for JJs. The first 
is the product IcRn, where Iq is the critical current 
and i?jv is the normal state resistance. IqRn should 
be of order A/e (where A is the superconducting gap 
of the leads and e is the charge of the electron) and in- 
dependent of device geometry^^. Two experiments^^'^^ 
have noted a smaller-than-expected value for the prod- 
uct IcRn in TI JJs. The second characteristic property 
is the "Fraunhofer-like" magnetic diffraction pattern, i.e. 
the decaying, oscillatory response of the supercurrent to 
the magnetic field B, applied perpendicular to the flow of 
the supercurrent. The first minimum in Iq should occur 
at B = Be, when one quantum of flux $o = h/2e (where 
h is Planck's constant) is passed through the area of the 
device, and recent reports on in TI JJs^'''^^ match this 
expectation. 

The main experimental results of this paper are a small 
value of IcRn that scales with 1/W, where W is the 
width of the JJ, and a value of Be that is ~5 times 
smaller than that expected from the device area. Nei- 
ther of these results is predicted or previously seen for 
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FIG. 1: Phenomenological model for Majorana fermions in topological insulator Josephson junction, a, Schematic of a 
topological insulator Josephson junction. Two superconducting leads are patterned on top of Bi2Se3 forming a junction with length L 
and width W. Along the width of the device, a one-dimensional wire of Majorana fermions results (purple), b, Energy levels near ip = it 
before momentum quantization along W (purple line) and after, where the topological state remains at E = (solid purple circle) and 
the first quantized energy level at the value Eq (empty purple circles), c. Current-phase relation resulting from momentum quantization, 
producing an anomalous peak at tp = n. d, Scanning-electron micrograph of a device similar to the ones measured in this Letter. 



conventional JJs nor TI JJs. To explain these experi- 
mental observations, we propose a two-fold phenomeno- 
logical extension to the model in Ref 4. First we take into 
account confinement along the ID Majorana wire, quan- 
tizing its energy levels at multiples of Eq = hvex/'^l, 
with i/ex the velocity of the carriers in the wire and I the 
length of the wire. In the present devices, the length of 
the wire is the width W of the JJ, hence Ec = hi/ex/^W. 
The effect of this quantization on the energy levels is 
shown in Fig. lb. If the E = state (Fig. lb, solid 
purple dot) is topological in nature, i.e. it is a neutral 
Majorana mode, such confinement should not affect its 
existence nor change its energy from zero. The contin- 
uum of energy levels E ^ 0, not protected from per- 
turbations, are quantized in multiples of Eq (Fig. lb, 
empty purple dots where only the first non-zero energy 
modes are shown for clarity). The second extension of 
Ref. 4 is to postulate that only the topological E = 
Majorana state carries supercurrent, while for E ^ 
(separated from the topological state by Eq) the trans- 
fer of carriers from one lead to the other is dissipative. 
This results in an IcRn set by the energy scale of con- 
fined modes along the width of the junction, rather than 
by A/e: Ic'Rn oc Ec/e = hVf,x/2eW . This also has 
an effect on the current-phase relation, which determines 
the supercurrent through the device as a function of Lp. 
The supercurrent-carrying E = Q state occurs at (/? = tt, 
producing a sharp peak in the current-phase relation at 
<p = TT (Fig. Ic), in contrast with the established sinu- 
soidal current-phase relation for conventional JJs^^, or a 
different-period sinusoid predicted for TI JJs^^. 



As described above, the result of confinement is to sep- 
arate in energy the £^ = ID modes (neutral Majorana 
modes) from the E ^ Q (charged) ID modes. A ID 
mode can carry a supercurrent, but the critical current is 
strongly suppressed by interactions between charges^*' 
Thus a neutral mode can carry a larger supercurrent than 
a charged mode. In our experiment the lowest-energy 
charged modes are accessed when the current through 
the neutral mode is 1 /iA. We estimate that the charged 
modes cannot carry this much supercurrent. Hence the 
charged modes in our devices act as resistors carrying 
current but not supercurrent. When the charged modes 
become energetically accessible (i.e. for energies > Ec) 
they shunt the junction, causing the JJ to undergo a su- 
perconducting to metallic phase transition^^, which takes 
places because the ID-charged-mode shunt resistance is 
greater than h/Ae^ . In this model, any additional su- 
percurrent through the bulk (which might be expected 
in existing TIs, given the finite bulk conductivity) also 
ceases when the shunt resistance of the surface become 
energetically accessible. 

To investigate the properties of JJs with TI weak links, 
Bi2Se3 crystals were synthesized either via a slow cooling 
of a binary melt of Bi and Se^^ or by vapor-liquid-solid 
(VLS) growth^^ and flakes of thicknesses between 50 and 
100 nm were obtained from these crystals via mechani- 
cal exfoliation^. Junctions of lengths L between 20 and 
80 nm and widths W between 0.5 and 3.2 /im were fab- 
ricated via electron-beam lithography, followed by depo- 
sition of Ti/AI electrical leads (see the Supp. Info for 
details of the sample fabrication). A scanning electron 
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FIG. 2: I-V response of a topological insulator Joseph- 
son junction. (main) V vs. I for a devices of dimensions 
(L, iy)={45nm, l/zm) for B=0, 2, 3, 5, 8, 10 mT and at a tem- 
perature of 12 mK. At B=0, Ic is 850 nA, wliicli is reduced upon 
increasing B. For this device, the product IqR^=30.6 fiV, much 
lower than theoretically expected for conventional JJs. (upper-left 
inset) I — V curves overlap for all values of B at V > 2A/e~300 /iV. 
(lower-right inset) Sweeps up (red) and down (black) in I show lit- 
tle hysteresis, indicating that the junction is in the overdamped 
regime. 



micrograph of a completed device is shown in Fig. Id. 
The DC response for a (L,W^) = (45 nm, 1 /im) junction 
at a temperature of 12 mK is shown in Fig. 2, where the 
DC voltage {V) is plotted as a function of the applied 
DC current (/). At i? = 0, a typical DC Josephson re- 
sponse is observed (red curve): for / < /c=850nA, V—0 
and a supercurrent flows. Applying B perpendicular to 
the top surface of the Bi2Se3 reduces Ic until 5=10 mT 
when the superconducting leads are driven normal and 
the I-V curve becomes linear. For I > Iq, there is an 
excess current due to Cooper pairs leaking into a low- 
barrier junction^"*; this excess decreases with B. For V > 
2A/e'-^300 /iV, all curves fall on top of each other (upper 
left inset of Fig. 2) for all values of i3. Absence of hystere- 
sis (lower right inset of Fig. 2), indicates that the junc- 
tion is overdamped, consistent with calculations (Supp 
Info). i?jv for this device is 35 and /ci?Ar=30.6 ^V. 
Theoreretical expectations^^ for diffusive and ballistic 
weak links predict Ic'Rn to be 281 /iV and 427 /iV, re- 
spectively, an order of magnitude higher than our mea- 
surements. As a control experiment, a device fabri- 
cated similarly to the TI JJs, except with a 75 nm-thick 
graphite weak link, has /(7i?Ar=244/iV (Supp. Info.), 
much closer to theoretical predictions. This suggests that 
something in the sample rather than the measurement 
setup reduces the values of Ic'Rn- 

Further insight into the nature of transport in TI JJs is 
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FIG. 3: 1/W dependence of IcRff. a, A comparison of two 
devices with similar Rf^ (56.1 and 51.5 O) and different widths 
W, 0.5 and 1 /im. The device with W=0.5 fj,m exhibits a larger 
Ic, in contrast to conventional JJs, where similar resistances lead 
to similar values of Iq- b, IcRn vs. 1/W for all 14 devices 
(synthesized via two methods: a binary melt and VLS) showing 
the general trend of IcRn 1/W. Using the phenomenological 
model of Fig. 1, a carrier velocity of 1.4it0.2x 10'* m/s is extracted. 



found by investigating the width dependence of the char- 
acteristic quantity IcRn- A comparison of two junctions 
with Rn = 56.1 and 51.5 O, and W=l and 0.5 fim, (both 
L=50 nm) is shown in Fig. 3a, demonstrating that halv- 
ing W doubles Ic despite similar values of Rn- Measure- 
ments of i?jv were carried out above the superconducting 
transition temperature of the leads in a four-terminal ge- 
ometry, eliminating the resistance of the cryostat lines, 
but not the contact resistance between Ti/Al and Bi2Se3, 
which varies from device to device without apparent cor- 
relation to geometry or effect on IcRn product. The 
mismatch of the I-V curves above Ic is due to the ex- 
cess current mentioned in connection with Fig. 2 (inset), 
which is typically of order /c^^ and hence larger in the 
W^=0.5 /xm device. The two curves approach each other 
as V approaches 2A/e (inset Fig. 3a). The values of 
IcRn for all 14 devices we measured that superconduct 
are shown as a function of 1/W in Fig. 3b. The trend is 
clear: a larger produces a smaller Jci? at. With benefit 
of hindsight, some previously reported experiments on TI 
JJs are consistent with IcRn being related to 1/W. In 
very wide junctions no supercurrent is observed at all^'^. 
For intermediate values of W similar to ours a low IcRn 
was measured^ ^. The velocity v^x of dissipative excita- 
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tions in the junction can be extracted from a linear fit of 
the data in Fig. 3b assuming a unity proportionahty con- 
stant between IcRn and Ec (i.e. IcRn=Ec /e), giving 
a value 1.4±0.2x 10"* m/s. Electrical transport measure- 
ments of the Fermi velocity vp for surface electrons in 
Bi2Se3 range from lO^m/s to lO^m/s bracketing 
the value i/ea;=4.2x 10^ m/s "^^ extracted from ARPES, 
and larger than our inferred value. In fact, a lower 
group velocity is expected for the bound pairs of electrons 
and holes that shuttle Cooper pairs across the device^^. 
Specifically for Majorana fermions - a subset of these 
bound pairs ~ the velocity vm has been predicted to be 
less than vp by a factor (A/^)^ in Ref. [4] and (A//i) 
in Ref. [28] (calculated for neutral modes created in 
graphene JJs), where ^ is the chemical potential of the 
TI weak link relative to the Dirac point of the surface 
states. The typical ratio for A//.t in our samples is 10"'^, 
giving an estimate for the Majorana velocity of 0.42 m/s 
for (A//.t)2 and 4.2x10^ m/s for (A/^). Given that we 
have ignored the proportionality constants in both the re- 
lationship between IcRn and Ec, and between vm and 
{A/fj,)vF, our results are reasonably near expectations 
for i^M given by Ref. [28]. 

The last characteristic response of TI JJs considered 
in this Letter is the magnetic diffraction pattern. The 
differential response dV/dI{B, I) for a {L,W)={55 nm, 
1.5 /um) device demonstrates an atypical relationship be- 
tween Ic and B (Fig. 4a). Three features are of note: Be 
is 5 times smaller than expected from the known device 
area; the shape of Ic{B) deviates from a typical Fraun- 
hofer pattern; and there exist added features in dV/dl 
at values of / above Ic whose magnitude and number 
depend on B. We consider the first two features below 
and address the third in dV/dl in the Supp. Info. The 
area of the devices is calculated as W * (L + 2\l), where 
Al=50 nm is the dirty London penetration depth for alu- 
minum (Supp. Info). The extracted Ic{B) is shown in 
Fig. 4b (red) and compared to the simulated Fraunhofcr 
pattern (black) for the device area (see Supp. Info, for 
information on the simulation). Be for this device is 
1.70 mT, whereas it should be 9.3 mT, based on the de- 
vice area measured from a scanning electron micrograph. 
We have measured a smaller-than-expected value of Be 
in all our devices. The three minima in Ic on each side 
of S = are unequally spaced, occurring at _B=1.70, 
6.25, and 11.80 mT. Even if the effective area of the junc- 
tion were larger for unknown reasons, fitting the central 
peak to a Fraunhofer pattern would produce minima at 
1.7, 3.4, and 5.1 mT, different from what is observed. 
The graphite control device exhibits a more conventional 
magnetic diffraction pattern (Supp. Info). 

Drastic narrowing of Ic (B) can occur if the current- 
phase relation (CPR) is peaked at certain values of ^p, 
as our phenomenological model suggests. We note that 
an anomalous, peaked CPR has been theoretically pre- 
dicted for a different device geometry, also a result of 
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FIG. 4: Magnetic interference pattern in topological in- 
sulator Josephson junction, a, Differential resistance dV/dl 
as a function of B and / showing an anomalous magnetic diffrac- 
tion pattern for a W=1.5 fim junction. Two features are of note: 
a smaller than expected value of Bq at 1.70 mT and a nonuni- 
form spacing between minima at values i?=1.70, 6.50, 11.80 mT. 
b, (main) Iq (B) (red) extracted from dV/dl in a is compared to 
the expected Fraunhofer pattern for the junction (black) where a 
reduction of the scale of the pattern and the nonuniform spacing 
are evident, (inset) A comparison of the simulated Frauhofer pat- 
tern for a sinusoidal (black) and an empirically-determined, peaked 
(red) CPR. The aperiodic minima observed in a are captured this 
CPR. 



the presence of Majorana fermions^^. In each case, one 
peak in the CPR occurs for each zero-energy crossing of 
states in the gap. For a single peak in the CPR, only 
a single, B=0 maximum is possible for |$| < $o (Supp. 
Info). The existence of multiple minima in Ic{B) for 
|"I>| < $0 indicates the presence of multiple peaks in 
the CPR - a result of coupling to fermionic modes in 
the device that create additional zero-energy crossings^°. 
For a given CPR, the resulting /c(B) is uniquely deter- 
mined; however the reverse is not true: a peak in the 
CPR at, for example, = 7r/2 produces the same Ic (B) 
as a peak at ip ^ Through simulation, we are able 
to show that a smoothly-varying CPR cannot capture 
our results (Supp. Info) and only a CPR with peaks can 
create magnetic diffraction patterns like those observed 
in our devices. The magnetic diffraction pattern for one 
such peaked CPR that also includes a conventional, si- 



5 



nusoidal contribution from the bulk is compared to the 
typical Fraunhofer pattern (black) in the inset of Fig. 4b. 
Some, but not all, of the features observed in the experi- 
ment arc captured by this CPR. Importantly it prodiiccs 
non- uniformly distributed minima at $o/4,3$o/4, and 
$0 (see Supp. Info for details) matching the aperiodic 
structure of the minima seen in the experiment. 
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Supplementary Information for Signatures of 

Majorana Fermions in Hybrid 
Superconductor- Topological Insulator Devices 

Materials and methods 

Single crystals, grown by the binary melt method, 
of (Bii_a;Sba;)2Se3 were synthesized by slow cooling 
a binary melt of Bi (99.9999%), Sb (99.999%) and 
Se (99.9999%) starting materials, mixed in the ratio 
0.33:0.05:0.62. The actual amount of Sb in the single 
crystals is approximately 0.01, as measured in micro- 
probe analysis. Bulk single crystals with dimensions 
IxlxO.lmm"^ showed a non-metallic resistivity at low 
temperatures with a bulk carrier density of 7x 10^^ cm~^ 
as deduced by Hall and Shubnikov-dc Haas analysis. 

Bi2Se3 nanowires, synthesized via the vapor-liquid- 
solid (VLS) method, were grown in a 12-inch horizontal 
tube furnace with a quartz tube. Bi2Sc3 source powder 
(99.999%, from Alfa Aesar) was placed in the center of 
the furnace (540 °C). The growth substrate, a Si wafer 
with a thermally-evaporated, 10 nm Au film, was placed 
at a downstream, lower temperature zone (350 °C). High 
purity Ar gas delivers vapor from the source materials to 
a growth substrate at 130 seem for 2 hours with ITorr 
pressure. 

Samples of 50-100 nm thickness were prepared by me- 
chanical exfoliation of the above material, using a method 
similar to exfoliation of graphene'^^. Prior to sample ex- 
foliation, SiO2(300nm)/Si wafers were cleaned in acetone 
and isopropanol. HMDS was deposited on the Si02 sur- 
face in efforts to reduce the amount of H2O contamina- 
tion. Samples were exfoliated using Nitto-Denko tape 
and immediately after exfoliation PMMA 950 A4 (Mi- 
crochem Crop.) was spun on the chip for 60 sees at 
4000 rpm, followed by a baked for 30mins at 180 °C. 
After electron beam exposure of the resist and devel- 
opment in MIBK/IPA 1:3, a 3min UV-ozone was used 
to remove residual PMMA in the exposure contact ar- 
eas. Ion milling of the exposed contact areas was per- 
formed to obtain low contact resistance to the Bi2Se3 
flakes. The recipe of the ion mill is a 300 V beam voltage, 
30 mA beam current producing an ion current density of 
0.1 mA/cm^, 60 V acceleration, 12 seem argon flow, and a 
10 sees exposure. Immediately following ion milling, the 
samples were loaded into c-beam evaporator with oil-free 
pumps and Ti/ Al (3 nm and either 60 nm or 100 nm) were 
evaporated at a rate of 1 A/s onto the sample at a pres- 
sure of ^5xl0~^ Torr at the start of each evaporation. 

Samples were measured in a ^He/^He dilution refrig- 
erator with a base temperature of 12 mK. To reduce the 
effects of thermal radiation, three stages of filtering were 
employed: room-temperature LC 7r-filters, a 3-pole RC 
filter anchored to the mixing chamber designed to filter 
in the frequency range of 1 MHz to 5GHz'^^, and lOkfi 
"on-chip" resistors used in conjunction with the inline 



capacitance of the twisted-pair wiring to produce a low- 
pass, RC-filter of cutoff frequency ~50KHz. All experi- 
ments were performed in an electrically- and acoustically- 
shielded room. Device DC resistance and differential re- 
sistance were measured using a standard DC and lock-in 
technique. The applied AC current was kept low enough 
such that the AC voltage across the device never ex- 
ceeded 500 nV, below the base temperature of the fridge 
k^T^l (UeV. To obtain accurate magnetic field values, a 
low-noise Yokogawa 7651 was used to supply the current 
to the superconducting magnet, producing a <0.3/iT 
noise level in the magnet. 

The transport properties of the exfoliated binary melt 
samples were characterized by making a device with con- 
tacts in a Hall-bar geometry and thickness of 100 nm. 
The extracted density, mobility and diffusion constant 
D of these samples are 7.8x 10^^ cm~^, 1950cm^/Vs, 
0.02m^/sec, respectively. We note that these values re- 
flect merged contributions from both the surface and 
bulk, where it is expected that the surface electrons have 
a much higher mobility'^^. Hence, the values for mobility 
and D are a lower bound for the values of the surface 
electrons. 

To extract the superconducting and normal state prop- 
erties of the leads, four samples were prepared with the 
same recipe as the leads of the J J devices: Ti/Al of 
thickness 3nm/(60, 100 nm), both of (L, I^)=(10,l) /xm. 
Each thickness of Al was placed either directly on the 
Si02 substrate or on top of a Bi2Se3 flake. Character- 
ization of these test samples was performed via a four- 
terminal resistance measurement at room temperature 
(290 K), 4K and base temperature of the dilution refrig- 
erator (12 mK). Little variation was found in the super- 
conducting properties of the Al patterned on the Si02 
substrate and on Bi2Se3. The resistivity of the samples 
at room temperature and 4K were 3.0x10"^ f2-m and 
0.3x 10~^ ri-m, respectively. The critical currents and 
magnetic fields for the 100 nm(60 nm)-thick Al-films were 
488M(208/iA) and 11(13) mT respectively. 

Characteristic lengths and regimes for JJs 

There are many different regimes in JJs that produce 
different behaviors. The regime applicable to the device 

is determined by the length of the junction L relative to 
five intrinsic length scales: coherence length of the Al 
leads ^Ai = l-6/im; the superconducting coherence length 
of the metallic (Bi2Sc3) weak link ^^r, and three prop- 
erties of the (Bi2Sc3) weak link unrelated to the super- 
conductor: the thermal length It^ the mean-free path 
le and the phase-coherence length (.^p. can be cal- 
culated from the diffusion constant by \/hD/A, where 
A is the superconducting gap of Al, which we calculate 
from the BCS equation A=1.76 fcBTc=151 /ieV, produc- 
ing a value for ^jv of ~ 280 nm. £t = ^/hD/2TTkBT, 
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which at 12 mK is 1.3 /zm. £e can be extract from D as 
D=l/2j/i?£e, giving ^e=80nm, where we have used the 
value of z^i?=4.2x 10^ m/s from ARPES measurements'^^. 

was estimated from the half-width at half-max of the 
weak antilocahzation correction to the longitudinal re- 
sistance (measured in the Hall bar) to be 650 nm. We 
note that £,n, ^t, and perhaps tip all represent lower 
bounds since we have used the average value of D from 
the surface and the bulk. The D for the surface elec- 
trons is expected to be higher. Further, W for the device 
should be compared to the Josephson penetration depth 
Aj — ^$o/(27r^ojc(2AL + L)), where Al is the London 
penetration depth ^ SOnm"^^. We place a lower bound 
on Aj by using the junction parameters that would make 
it the smallest, i.e. using the largest measured value of 
jc and the longest value of L: producing a value for 
ynin q£ IQ /im. This puts our junctions in the regime 
L^^e < £,N,£.Ai,i^ip,LT and W < Xj, i.e. the short 
junction, quasi-ballistic regime. In this regime, IcRn 
is expected to be described by either the KO-1 or KO-2 
theory^^ and independent of W, neither of which expec- 
tation matches our data. 



Supercurrent carried by the bulk 

The typical thickness t for our devices is 75 nm. The 
calculated ^at is larger than this value, which means that 
both the surface and entire bulk should experience the 
proximity effect. We expect, however, that the supercur- 
rent carried by the bulk to be smaller than that of the 
surface for two reasons. First, the Cooper pairs have to 
travel a longer distance from left to right lead through 
the bulk than through the surface. We can estimate this 
additional length by assuming the Cooper pairs have to 
go on average a distance t/2 from the left lead to the 
bulk, then another t/2 going from the bulk to the right 
lead. This makes the effective length of the device for the 
bulk states to be L^uik = L + t. The second is the lower 
mobility of the bulk, producing a smaller value of D. For 
the device in Fig. 4 of the main manuscript L—bb nm 
and t=95nm, making Lt,„;fe=150 nm. To calculate the 
reduction of the critical current, we assume D for the 
bulk states is 10 times smaller than the calculated value 
(a factor 12 has been experimentally determined'^^ for 
the reduction of D for the bulk), giving a A^"^' ~ 65 nm. 
This estimate yields a bulk critical current of order five 
times smaller than that of the surface state"^^. 
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FIG. 5; Additional features in dV/dl as function of temper- 
ature, Idc I and B. a, Differential resistance dV/dl as a function 
of temperatures Tf as Tf is lower through Tq of the leads. For this 
device, a sharp transition to the superconducting state is observed 
at Tj-=850mK. b, Additional features in dV/dl as a function of / 
are observed above Ic- c, dV/dl for values of Tf=12, 75, 150, 300, 
400 mK. Shoulder dips are observed in dV/dl for the two temper- 
atures 300 and 400 mK. d, dV/dl extracted for B=0, 0.75, 2.50, 
8.00 mT from Fig. 4a of the main text. Shoulder dips are also 
apparent when the magnetic field is increased to 0.75 mT. 



type of junction is determined the quality factor Q = 
y/2eIcC/h Rn , where C is the capacitance of the J J. We 
estimate C to be 0.5pF using a parallel plate capacitor 
model between the entire area of the leads of the device 
and the degenerately-doped Si, yielding (5=1.2. This is 
neither overdamped or underdamped, but at an inter- 
mediate value. Numerical evaluation of the /- V charac- 
teristics in this regime shows results similar to those of 
an overdamped junction The effect of thermal radia- 
tion on overdamped junctions has been calculated'^*, and 
shown to cause a smearing of the transition between the 
superconducting and normal state. Since the transitions 
we measure are still very sharp, we can rule out thermal 
fluctuations being a source of the reduction of IcRn- We 
can also calculate what the effect of thermal fluctuations 
would be had the junctions been in the underdamped 
regime the effective thermal radiation temperature 
would have to be 3.4 K to reduce IcRn by the factors 
we observe. This high temperature is unphysical in our 
well-filtered setup. 



Mechanisms for the reduction of IcRn 

We next consider the other possible mechanisms for 
the reduction of IcRn- A typical reduction comes from 
thermal noise, where the effect on IcRn depends on the 
type of junction (overdamped or underdamped). The 



Additional features in the Differential Response 
dV/dl as a function of T and B 

The differential response dV/dl for a (L, W^)=(50 nm, 
0.9 /im) device as a function of fridge temperature Tf 
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is shown in Fig. Sla, where the device exhibits a 
sharp transition to dV/dI=0 below Ty=850mK of the 
leads. Above Tc, the devices shows a metallic tem- 
perature dependence of resistance. The differential re- 
sponse dV/dI{I) of this junction, typical for all devices 
measured, at r/=12mK is shown in Fig. Sib. Even 
aside from the transition from a superconducting state 
to a normal state, which produce large peaks at / = Jc, 
several features are evident and are indicated by labels 
1-4. The first occurs in all devices as either a peak 
or a sharp rise in dV/dl (like the one shown in Fig. 
Sib), at V — ^\ . The second, broad peak appears at 
V — 300 /^V, near 2A/e of Al, where a peak in resistance 
is expected due to the quasiparticle contribution to the 
resistance from the leads'^^. The third and fourth peaks 
are both features occurring above 2A of the leads, occur- 
ring at = 500 and V — 655 /^V, respectively. Peaks 
above 2 A have also been seen previously in TI JJs'^^". 

The effect of Tj on dV/dl is shown in Fig. Sic. A 
typical trend is observed, where increasing Tf reduces 
Ic- Additional "shoulder" dips appearing in dV/dl are 
evident in the T/=300mK and 400 mK traces. On the 
400 mK trace, two dips are seen, at 7=500 nA and 620 nA 
respectively. 

Additional shoulder features are also seen in dV/dl 
above Ic when a magnetic field B is applied (Fig. Sid). 
These are revealed in cuts of constant B in dVdI{B, I) 
for the device shown in Fig. 4 of the main text. At B = 0, 
there are two peaks in dV/dl at I = ±/c, and no other 
apparent features. At _B=0.75mT, dips are evident at 
/ = ±1.6 /xA, T^=12.5/iV above the resistive transition. 
The single dip occurs for 0.6<-B<l mT. Two dips appear 
beginning at B^l mT and persist throughout the entire 
second and third lobes, although to a much weaker degree 
in the third lobe. The double-dip is evident in the cut 
at i?=2.5mT, occurring at /±=0.9 and 1.2 /iA and a 
value of V=7A and 12.1 /iV. The dips at 5=8 mT have 
all but disappeared, remaining weakly at /±=0.5 and 
0.9 /iA and V—6.3 and 12.1 /iV. Currently, we have no 
explanation for these features which appear in all devices 
measured. 



Andreev Bound States and the Current-Phase 
Relation 

In this section we consider qualitatively the Andreev 
bound states (ABS) and the resulting current-phase re- 
lation (CPR). States in a J J within the superconducting 
gap are associated with the coupled electron-hole pairs 
that are responsible for transferring Cooper pairs from 
one lead to the other'^-'^^. These are known as ABS. To 
describe the CPR in our junction, we consider several 
different energy spectrum of ABS. The energy levels of 
conventional ABS are E{(p) cx ±Ay/l -Tsin^{(p/2) 
where T is the transmission probability of Cooper pairs 
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FIG. 6: Determining the Current-Phase Relation from the 
Andreev Bound State Energy Spectrum. The CPR can be 

determined from the Andreev bound state (ABS) energy spectrum 
by taking the derivative of E with respect to ip. Qualitatively 
shown is the ABS spectrum (a) and the corresponding CPR (b) for 
conventional, ballistic JJs with a transmission probability T=0.9. 
The same plots (c,d) are shown for a, p + ip and TI (without the 
lateral confinement of our model) JJs. ABS spectrum and CPR 
(e,f ) for the phenomenological model of the main text. The inset 
of e shows the dependence oiT on(p in our model. The CPR of (f) 
is peaked at the values of ip where the zero-energy crossings occur, 
g, The width of T in the inset of e and the corresponding CPR are 
set by Ec , the lateral confinement energy in the phenomenological 
model, h, (inset) ABS spectrum similar to S2e except with a anti- 
crossing in the spectrum at ip = n. (main) CPR resulting from the 
ABS of the inset. 



from left to right lead. From the energy spectrum of the 
ABS, the current phase relation (CPR) can be calculated 
as I{(p) (x dE/dcp^^'^ . The energy spectrum for the ABS 
in a conventional JJ with T=0.9 is shown in Fig. S2a 
for the occupied (solid) and unoccupied (dashed) ABS 
states. The corresponding CPR, calculated by differen- 
tiating the E{(f) plot, is shown in Fig. S2b, where a 
saw-tooth behavior is seen, occurring with a period of 
27r. Note that for low values of T, this CPR becomes 
the conventional sin {ipf^\ For p + ip and TI JJs, the 
energy of the ABS becomes E{(p) (x ztAcos{(p/2), pro- 
ducing a zero-energy crossing at iy9 = tt'^^^. There is a 
correspondence between the p + ip and conventional en- 
ergy spectrum of ABS for T = 1. The corresponding 
CPR is I{ip) (X sin{tp/2), resulting in an anomalous Att- 
periodic CPR, predicted for p + ip JJs'^^^ and TI'^^^ JJs. 
The energy spectrum of the ABS and the CPR in this 
case are shown in Fig. S2c,d. 
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In our phenomenological model, the transmission co- 
efficient T depends on Lp, where the ID charge modes at 
if ^ Ti have a low value for T - a result of the inter- 
actions in ID'^^* - and the neutral modes have a much 
higher value for T. This T{ip) produces a CPR that is 
peaked at values of tf where the zero-energy modes oc- 
cur. Fig. S2e shows an ABS spectrum for a T{ip) (inset 
of Fig. S2e), corresponding to T{ip)=l for ip = it and 
T(lp)—1/1Q for (/5 7^ TT. The exact value for T away from 
if = TT can be changed without much change in the re- 
sults . To connect the high- and low-T regions, we assume 
that the width of the peak in T{ip) near tt is set by the 
energy scale £'c, the energy associated with momentum 
quantization in the phenomenological model. This sce- 
nario is shown in Fig. S2g, where the linearized E{ip) 
is calculated from Eq. 4 of Ref. [SI 5]. Using a value 
of Ec—bIcRn—'^^ A^eV for the device in Fig. 4, values 
of 167° and 193° for (fi and Lp2 are obtained. T then 
is modeled as a gaussian with a half max at the values 
ipi and ip2- The calculated ABS spectrum for this trans- 
mission coefficient is shown in Fig. S2e, where a tran- 
sition through E = occurs on the scale of ipi and ip2 
around ip = tt. The corresponding I{(p) is shown in Fig. 
S2f, where is it seen that the iy9-dependent T produces a 
anomalous CPR that is peaked around the values of (p 
where the zero-energy modes occur. This CPR is also 
Att periodic. A similar CPR was obtained in Ref. S16, 
where a peak in the CPR was also found to be correlated 
with the zero-energy crossing in the ABS. In Fig. S2c 
and S2e, it was assumed that sermonic parity was con- 
served, which results in a protected crossing of the two 
states (solid and dashed lines) at (p = tt"^^^. If the parity 
of the junction is not preserved, an anti-crosssing of these 
two states occurs (inset of Fig. S2h), producing a CPR 
shown in Fig. S2h. 



Simulation of Josephson Effect in the Presence of 
Magnetic Flux 

The critical current through the devices as a function 
of applied field was calculated in a manner closely fol- 
lowing that of Tinkham'^^, but with the extended junc- 
tion model modified to allow non-standard current-phase 
relationships and spatial inhomogeneities critical current 
density. With the flake surface parallel to the x — y plane, 
applied field B along z, and current along cc, the phase 
difference between the leads as a function of the position 
y along the junction is given by 



2tt 

$0 



y rL+2XL 

dy / dxB^{x,y), 
Jo 



(1) 



where ipg is the phase difference at y=0, $o is the mag- 
netic field flux quantum, L is the length of the junction, 
and Xl is the London penetration depth. This follows 
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FIG. 7: Magnetic diffraction patterns for conventional 
current-phase relations. Magnetic diffraction patterns (red) for 
a, sinusoidal /(</p), b saw-toothed sinusoidal I{<f), c,d, sinusoidal 
I{(p) with two different nonuniform critical current line densities 
Kc- In b-d, the magnetic diffraction pattern for the sinusoidal 
I{ip) (black) from a is plotted for comparison. In each figure, the 
current-phase relations used to produced the magnetic diffraction 
pattern is shown in the inset. 



from integrating the vector potential A around a rect- 
angular contour that includes a path at y'=0 and y'=y 
and taking advantage of the fact that A ~ (<I>o/27r) V7 in- 
side the superconducting leads, where 7 is the non-gauge- 
invariant phase [defined so that (p = A7 — (27r/<I>o) / A • 
ds]. With normalized current-phase relationship (CPR) 
i{(p) = I{ip)/Ic and critical current line density Kc{y), 
we obtain critical current 



Ir 



W 



dyKciy)i{(p{y)). 



(2) 



Note that this analysis is equally valid for 27r- and 47r- 
periodic phase relations, so long as for the latter case we 
allow Lp{y) to range up to 47r and define i{(p) over the full 
range. For a typical 27r-periodic current-phase relation- 
ship I{(p) = Icsin{(p) and uniform current density and 
field, the critical current drops to zero when an integer 
number of flux quanta are threaded through the junc- 
tion area, corresponding to an integer number of cycles 
of Josephson current. A simulation of magnetic diffrac- 
tion pattern (MDP) for i{ip) ~ sin{(p) shows the typical 
Fraunhofer pattern'^^ (Fig- S3a). Deviations from this 
typical pattern can occur because of non-sinusoidal CPR 
or from nonuniform critical current distribution along the 
junction. For a ballistic JJ, the CPR becomes more of a 
saw-toothed shape, for which the corresponding MDP is 
shown in Fig. S3b. There is a small change in the shape 
of the MDP, but the minima still occur at integer multi- 
ples of <I>o. Aside from a changing CPR, deviations from 
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FIG. 8: Magnetic diffraction patterns for anomalous 
current-phase relations. Magnetic diffraction patterns (red) 
for a two gaussian peaks in at (/p = itTr. b, tlie same as a 

but witii an added conventional sinusoidal of 1/5 the am- 

plitude of the anomalous peaks, c, the same a, with additional 
peaks at (/p = 7r/2, 37r/2 and 57r/2,77r/2, resulting from additional 
zero-energy crossings, d, the same as c but with a conventional 
sinusoidal /(v) of 1/5 the amplitude of the anomalous peaks. In 
each figure, the current-phase relations used to produced the mag- 
netic diffraction pattern are shown in the inset and each diffraction 
pattern is compared to the conventional sinusoidal case (black), e, 
A peaked CPR for an ABS spectrum in Fig. S2h, with additional 
zero-energy crossings at i/p = 7r/2, tt and 37r/2. f, same as e except 
with the sinusoidal bulk contribution added. 

the Fraunhofer pattern for short junctions can occur for 
a non-uniform current distribution. Reductions of the 
the location of the first minimum as large as a factor of 2 
can occur in the extreme limit when all of the current is 
concentrated near the end of the junction"^^^. The MDP 
in this extreme case is shown for two different current 
distributions in Fig. S3c, d. None of these deviations 
produce a MDP similar to Fig. 4 of the main text. 

Using this simulation, we empirically determine a CPR 
to match the MDP we observe in our devices. Not 
any CPR is possible: the CPR must be antisymmet- 
ric I{'~p) = —I (—if), it must come from the deriva- 
tive of the energy of the junction and must average to 
zero over a period such that the energy of the junction 
E{ip) oc J dipl{ip) does not continue to grow as a func- 



tion of (p . In addition, the CPR must contain a term 
for the contribution of the bulk, which was estimated 
in our devices to be smaller than the surface, but not 
zero. In the following, we assume the bulk has a CPR 
I{ip) (X sin{(f). To begin, we assume the simplest CPR 
that is peaked at values of tt and Stt, with alternating 
signs but equal amplitude to satisfy the conditions on 
the CPR. This result is shown in Fig. S4a, where we 
see that the inclusion of peaks at ip — tt, Stt produces 
a narrowed MDP. It does not, however, reproduce the 
additional minima seen in Fig. 4. Note that this CPR 
is 47r-periodic, as expected for TI JJs. The first correc- 
tion to the simple, peaked CPR above comes from the 
addition of the bulk term, which we include as a simple 
sinusoidal dependence on ip, sin{ip), with 1/5 the ampli- 
tude (as estimated above) of the peaks in the CPR and is 
27r-periodic. Compared to Fig. S4a, additional features 
in the MDP, are thereby introduced at n^o, where n is 
an integer. To get closer to the observed MDP, we add 
more zero-energy crossings, predicted in Ref. [S17]. We 
add two additional zero-energy modes or either side of 
ip = TT and Stt, at 1^9 = 7r/2, 37r/2 and 5n/2, 77r/2, and the 
corresponding MDP is shown in Fig. S4c,d. The MDP is 
calculated without (S4c) and with (S4d) the sin{(p) con- 
tribution from the bulk. It is seen in Fig. S4d, a small 
notch in the MDP appears below the first expected min- 
ima, occurring at a value of 3$o- A CPR similar to that 
used in Fig. S3c,d was calculated as an anomalous contri- 
bution to the CPR from Majorana fermions'^^^, a result 
of difference in Josephson current from the even and odd 
Majorana occupancies in the JJ. The last CPR consid- 
ered are of the form from Fig. S2h and are shown in Fig. 
S4e,f, where we again have used the multiple zero-energy 
crossings of Ref. [SI 7]. The bulk importance of the in- 
clusion of a bulk term is evident here, where without the 
sinusoidal addition no sub-$o dips occur (Fig. S4e). The 
MDP of Fig. S4f produces minima at values $o/4, 3<I'o/4 
and $0, near the experimentally observed minima. We 
note that we cannot determine a unique CPR from the 
diffraction pattern but can only supply evidence for a 
peaked CPR. In our simulations, we were unable to ob- 
tain sub-$o dips for any uniform (i.e. non-peaked) CPR, 
but were only able to obtain MDPs that look similar to 
our experimental results using peaked CPRs. 

MDPs for a varying bulk contribution are shown in 
Fig. S5. Fig. S5a is a 2D plot of the MDP for the CPR 
of Fig. S4e as a function of $ and the amplitude of the 
sinusoidal contribution (Sine Amp). When Sine Amp=0 
(all surface), the resulting MDP is the same as Fig. S4e, 
while for Sine Amp=l (all bulk) the MDP is same as Fig. 
S3a. Fig. S5a shows that sub-$o minima are obtained for 
values of Sine Amp less than ~0.4. For Sine Amp=0.5 
(red), the MDP obtained is shown in Fig. S5b, where 
it is seen that a narrowing of the central feature occurs 
when compared to a Fraunhofer pattern (black) , but the 
minima still occur at integer values of $0. A MDP similar 
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FIG. 9: Dependence of the MDP on the bulk contribution. 

(a) A 2D plot of the MDP as a function of $ and the amplitude of 
the sinusoidal bulk contribution (Sine Amp). Sub-<i>o minima are 
observed for values of Sine Amp less than ~0.4. (b) A ID cut in a 
at a value of Sine Amp=0.5 (red), showing a narrowed central lobe 
[compared to the Fraunhofer pattern (black)] in the MDP but no 
minima below <l>o. 



to the one calculated in Fig. S5b is seen in Fig. 2b of 
Ref. [S20]. 



A Graphite Control Device 

To confirm that the anomalous features observed in 
TI JJs are a result of the interface between a TI and 
a conventional superconductor, a JJ was created using 
75 nm-thick graphite as the weak link. The details of the 
device fabrication are identical to the fabrication of TI 
JJs, except the Ar ion mill step was not needed to pro- 
duce low contact resistance. Instead, a 5min UV-ozone 
exposure, typical in the creation of graphene devices'^^^, 
was used. The results, taken at r/=12mK, are shown 
in Fig. S5 where it is seen that the value of IcRn is 
much higher at 244 /^V (Fig. S5a) and Be is -13 mT, 
closer to the predicted value of 14.3 mT (black line) for a 
device of {L, iy)=(0.05,l) iim. The control device allows 
for the exclusion of two factors that might have caused 
the anomalous results in TI JJs. In Fig. S5a, a plot of V 
vs. / shows a value for the critical current of 8.1 /iA on 
a device of resistance 30.2 fi. The larger value of IcRn 
obtained for the control device rules out the possibility 
of the thermal effects caused by poor filtering being a 
source of the reduced IcRn in the TI JJs. In Fig. S5b, 
Ic vs. B (red), extracted from a plot of dV/dI{B, T) (in- 
set of S5b) shows a more conventional dependence on B, 
closely matching the expected pattern (black). The high 
value of Be in the control device rules out the possibility 
of flux focusing'^^^ reducing the values of Be in TI JJs. 
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FIG. 10: Superconducting properties of a graphite control 
sample, (a) Vvs. / for the (L,W)={50nm,l fim) control device 
created from a 75 nm-thick graphite weak link, producing a IcRn 
of 244 /iV, larger than the observed in TI JJs. (b) Extracted Iq 
vs. B (red) from dV/dl (inset) giving a value of Be of 13 mT, close 
to the value of 14.3 mT expected from the device area. 
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